 &MODEL
 NLIST   =           1,
 NSTEP   =        2000,
 EPSI    =  1.0000000E-07,
 WGHT    =   1.000000    ,
 NOBS    =         120,
 NVAXZ   =           4,
 NVAX    =           4,
 NORDER  =           2,
 NVAZ    =           0,
 NLAMZ   =           0,
 NLAMY   =           1,
 NLAMX   =           1,
 NSCALE  =           0,
 NEXPY   =           0,
 NRES    =           0,
 NCOR    =           0,
 NCOV    =           0,
 VARNAME = LOG+A1+A2,
 NVARIANT        =           5
 /
 &AUTOCOR
 RHO     = 2*0.000000000000000E+000  ,
 PI      = 2*1.00000000000000        ,
 NTPI    = 2*0,
 MATFILE = R1.in                                                                           R2.in                                     
                                       
 /
 &YBOXCOX
 VLY     =  0.000000000000000E+000,
 NTLY    =           0,
 NUMLXY  =           1
 /
 &XBOXCOX
 NUMX    =           1,           2,           3,           4, 46*0,
 NUMLX   = 4*1, 46*0,
 VLX     = 10*0.000000000000000E+000  ,
 NTLX    = 10*0
 /
 &YXZDATA
 YXZFILE = Canada1976.in                                                                   ,
 YXZFMT  = (5e15.5)                                                                        
 /

VARIABLE       MINIMUM         MEAN          MAXIMUM   NULL OBS.      SCALE
DEPENDENT
TOTFLOW     0.2577000D+04  0.1476588D+06  0.2736000D+07     0     0.1000000D-05
INDEPENDENT
POP         0.8591300D+05  0.5297039D+06  0.2802800D+07     0     0.1000000D-05
REVE        0.3698300D+04  0.4258776D+04  0.4808200D+04     0     0.1000000D-02
LANGFR      0.5600000D+01  0.7329333D+02  0.9990000D+02     0     0.1000000D+00
UTILITY     0.5655600D-04  0.4915871D-01  0.1172800D+01     0     0.1000000D+01
CONSTANT    0.1000000D+01  0.1000000D+01  0.1000000D+01     0     0.1000000D+01


MRS of X(l) with respect to X(k) - at the sample means
===================================================================================================================================
        X(k) POP         REVE        LANGFR      UTILITY     CONSTANT  
X(l)
===================================================================================================================================

POP          0.1000D+01  0.7451D+02  0.3456D+03  0.4914D+07 -0.4112D+07

REVE         0.1342D-01  0.1000D+01  0.4639D+01  0.6595D+05 -0.5518D+05

LANGFR       0.2893D-02  0.2156D+00  0.1000D+01  0.1422D+05 -0.1190D+05

UTILITY      0.2035D-06  0.1516D-04  0.7034D-04  0.1000D+01 -0.8368D+00

CONSTANT    -0.2432D-06 -0.1812D-04 -0.8406D-04 -0.1195D+01  0.1000D+01


INITIAL ESTIMATES

             LAM    PARAMETER      STD-ERROR        STUDENT-T    GRADIENT       DERIV.Y         ELAST.Y
                                                                                DERIV.E(Y)      ELAST.E(Y)


POP           X1  0.1389205D+01  0.5457770D-01     25.4537107 -0.1024327D-05  0.3872510D+00      1.3892049
                                                                              0.4970914D+00      1.3892049

REVE          X1  0.8322593D+00  0.6932686D+00      1.2004861 -0.6618732D-06  0.2885581D+02      0.8322593
                                                                              0.3626500D+02      0.8322593

LANGFR        X1  0.6643968D-01  0.5550501D-01      1.1970033 -0.3276004D-06  0.1338513D+03      0.0664397
                                                                              0.1916400D+03      0.0664397

UTILITY       X1  0.6335200D+00  0.1496708D-01     42.3275532  0.4788431D-06  0.1902915D+07      0.6335200
                                                                              0.2843971D+08      0.6335200

CONSTANT      X0 -0.1078358D+02  0.5515421D+01     -1.9551693 -0.7920692D-07 -0.1592291D+07    -10.7835828
                                                                             -0.2006817D+07    -10.7835828

VAR(W)            0.1312324D+00  0.1694202D-01      7.7459667  0.0000000D+00

INITIAL VALUES
 
RHO1           0.0000000D+00
RHO2           0.0000000D+00
PI1            0.1000000D+01
PI2            0.1000000D+01
 
LAM(Y)         0.0000000D+00
LAM(X1)        0.0000000D+00


LOG-LIKELIHOOD = -0.1318235D+04

NUMBER OF OBSERVATIONS              =            120

PEARSON-R2(E) (unadjusted)          =         0.8850
PEARSON-R2(E) (adjusted for D.F.)   =         0.8810

PSEUDO-R2(L)  (unadjusted)          =         0.9441
PSEUDO-R2(L)  (adjusted for D.F.)   =         0.9422

PROBABILITY of Y to be at the limit =         0.0000
MEAN VALUE of E(Y)                  =  0.1860993D+06

MEAN of residual W                  = -0.8662091D-10
VARIANCE of residual W              =  0.1312324D+00
95% CONFIDENCE INTERVAL for residual W
Lower Bound = -0.7100297D+00   Upper Bound =  0.7100297D+00




initialization cards for control routine follow - -  id card - no. parameters , no. steps , step size 
minimization tolerance - starting values of parameters - relative uncertainties of parameters (error estimates) -  
and starting direction  (relative parameter increments for 1st step)
identification line
         2      2000   0.1E+01   0.1E-06
  0.00000E+00  0.00000E+00
  0.10000E+01  0.10000E+01
  0.00000E+00  0.00000E+00

first entry to fcn , parameters of current step follow - - f,x(i)
   0.1318235024D+04  0.00000D+00  0.00000D+00

new absolute minimum , parameters for fcn entry     6  follow - - f,x(i)
   0.1306134110D+04  0.43801D+00  0.34447D+00

new absolute minimum , parameters for fcn entry     8  follow - - f,x(i)
   0.1306119357D+04  0.45287D+00  0.32684D+00

new absolute minimum , parameters for fcn entry    10  follow - - f,x(i)
   0.1306119328D+04  0.45343D+00  0.32722D+00

new absolute minimum , parameters for fcn entry    11  follow - - f,x(i)
   0.1306119328D+04  0.45344D+00  0.32722D+00

new absolute minimum , parameters for fcn entry    12  follow - - f,x(i)
   0.1306119328D+04  0.45344D+00  0.32722D+00

new absolute minimum , parameters for fcn entry    13  follow - - f,x(i)
   0.1306119328D+04  0.45344D+00  0.32722D+00


the function (fcn) has not improved by    0.10000D-06 after 3 steps which include a minimum - - -


therefore the run has been terminated
   0.1306119328D+04  0.45344D+00  0.32722D+00


all numbers of last entry refer to final pt. of run


the lowest value of the function (fcn) was   0.1306119328D+04 for entry   13 - - x(i) follows
  0.45344E+00  0.32722E+00



run statistics follow - -

 nfcn=   13 nst=    6


MRS of X(l) with respect to X(k) - at the sample means
===================================================================================================================================
        X(k) POP         REVE        LANGFR      UTILITY     CONSTANT  
X(l)
===================================================================================================================================

POP          0.1000D+01  0.1162D+03  0.1222D+04  0.4542D+07 -0.6068D+07

REVE         0.8609D-02  0.1000D+01  0.1052D+02  0.3911D+05 -0.5224D+05

LANGFR       0.8181D-03  0.9503D-01  0.1000D+01  0.3716D+04 -0.4964D+04

UTILITY      0.2201D-06  0.2557D-04  0.2691D-03  0.1000D+01 -0.1336D+01

CONSTANT    -0.1648D-06 -0.1914D-04 -0.2014D-03 -0.7486D+00  0.1000D+01


FINAL ESTIMATES

             LAM    PARAMETER      STD-ERROR        STUDENT-T      STUDENT-T    GRADIENT       DERIV.Y         ELAST.Y
                                   (CONDIT.)        (CONDIT.)      (CONDIT.)                   DERIV.E(Y)      ELAST.E(Y)
                                                    (PAR = 0)      (PAR = 1)

POP           X1  0.1413609D+01  0.5518806D-01     25.6144048                 0.1509000D-07  0.3940539D+00      1.4136093
                                 0.5518806D-01     25.6144048                                0.4525989D+00      1.4136093

REVE          X1  0.1320178D+01  0.8443306D+00      1.5635796                 0.9803065D-08  0.4577276D+02      1.3201781
                                 0.8443306D+00      1.5635796                                0.5289322D+02      1.3201781

LANGFR        X1  0.2390841D+00  0.9989701D-01      2.3933058                 0.4985512D-08  0.4816655D+03      0.2390841
                                 0.9989701D-01      2.3933058                                0.6463981D+03      0.2390841

UTILITY       X1  0.5959098D+00  0.2803633D-01     21.2549097                -0.8578515D-08  0.1789944D+07      0.5959098
                                 0.2803633D-01     21.2549097                                0.2657740D+08      0.5959098

CONSTANT      X0 -0.1619336D+02  0.7041822D+01     -2.2995979                 0.1169013D-08 -0.2391092D+07    -16.1933597
                                 0.7041822D+01     -2.2995979                               -0.2775317D+07    -16.1933597

VAR(W)            0.9495319D-01  0.1295236D-01      7.3309562                 0.0000000D+00
                                 0.1295236D-01      7.3309562
 
RHO1      (E)     0.4247198D+00  0.1255040D+00      3.3841145                 0.1133974D-04
                                 0.1255040D+00      3.3841145
RHO2      (E)     0.3160229D+00  0.1685477D+00      1.8749765                -0.1647057D-04
                                 0.1685477D+00      1.8749765
PI1       (F)     0.1000000D+01
PI2       (F)     0.1000000D+01
 
LAM(Y)    (F)     0.0000000D+00
LAM(X1)   (F)     0.0000000D+00


LOG-LIKELIHOOD = -0.1306119D+04

NUMBER OF OBSERVATIONS              =            120

PEARSON-R2(E) (unadjusted)          =         0.9202
PEARSON-R2(E) (adjusted for D.F.)   =         0.9159

PSEUDO-R2(L)  (unadjusted)          =         0.9771
PSEUDO-R2(L)  (adjusted for D.F.)   =         0.9759

PROBABILITY of Y to be at the limit =         0.0000
MEAN VALUE of E(Y)                  =  0.1713861D+06

MEAN of residual W                  =  0.4623180D-02
VARIANCE of residual W              =  0.9495319D-01
95% CONFIDENCE INTERVAL for residual W
Lower Bound = -0.5993405D+00   Upper Bound =  0.6085869D+00


CPU TIME:      0.25 sec
